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This paper describes an approach to construct temporally shaped control pulses that drive a 
quantum system towards desired properties. A parametrization in terms of periodic functions with 
pre-defined frequencies results in a smooth, typically simple shape of the pulses. This enables 
implementation of control without broadband sources and holds the potential to improve our un- 
derstanding of the functionality of optimal control. 
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Research on quantum mechanical systems is currently 
undergoing a process of substantial changes: whereas in 
the last decades the effort was mostly on the observa- 
tion and description of quantum mechanical properties, 
currently the option to manipulate and control quantum 
systems is moving into focus. On the one hand, this is 
due to the technological advances that permit isolating 
quantum systems e.g. in ion traps [1] [2] or optical lat- 
tices [3] , manipulating them coherently and letting quan- 
tum systems of different kinds interact with each other 
[H[S]. On the other hand, there is the prospect to exploit 
intrinsically quantum mechanical properties to engineer 
devices with performance characteristics far beyond the 
classically achievable. Whereas secure communication [B] 
or teleportation [7j based on quantum protocols is well 
established by now and the possibility to simulate quan- 
tum mechanical many-body systems with quantum simu- 
lators [SJIH] is a growing field, new perspectives to exploit 
quantum mechanical coherence phenomena, for example 
in energy provision [TOlIII] are just emerging. 

Beyond the experimental capability to control quan- 
tum systems, any type of quantum engineering also needs 
appropriate theoretical tools that allow an experimental- 
ist to extract optimal performance under given limita- 
tions of control, as typically imposed by power or fre- 
quency range of driving fields. Optimal control theory 
[T2TfT5] provides elaborate and efficient schemes to iden- 
tify e.g. shapes of laser or mircrowave pulses that drive 
a system towards desired properties. A particularly as- 
tonishing property of pulses designed by optimal control 
techniques is their robustness against experimental im- 
perfections [16»11Tj : for example inhomogeneous broaden- 
ing in ensembles of quantum systems can be compensated 
essentially completely through suitably designed control 

pulses PSH23]- 

A disadvantage of these pulses is that they typically 
contain many frequency components; besides potential 
experimental challenges to generate such pulses, the com- 
plicated structure of these pulses render it essentially im- 
possible to understand why they result in their astonish- 
ing performance. In particular if we want to push the 
envelope to large many-body systems, the answer to the 
question of l why l will become more and more important 
rather than the observation 'that' one can identify suit- 
able pulses. The aim here is therefore to strive for an 



approach that permits to limit a pulse to given frequency 
components [31], resulting in rather simple pulse shapes. 

The typical scenario of optimal control theory is dy- 
namics induced by a system- or drift-Hamiltonian T-Ld 
and an often time-dependent control Hamiltonian Hc(t). 
The system Hamiltonian is assumed to be inaccessible 
to control, i.e. it contains no adjustable components, 
whereas the control Hamiltonian can be tailored in a 
time-dependent fashion. The central question of inter- 
est is to determine how the propagator U induced by the 
time-dependent Hamltionian % = %d + He (t) changes 
if the control Hamiltonian is varied, since this provides 
the basis for an iterative improvement of He (t) [H E] • 
Typically, the time-dependence of the control Hamilto- 
nian stems from the time-dependence of an externally ap- 
plied control field, like a laser- or microwave field, so that 
Hc{t) = fi(t)'Hi, where the time-independent Hamil- 
tonians Hi describe the coupling of the system to the 
various control fields, and fi(t) is the time-dependence 
of the amplitude of the corresponding field. The objects 
of interest thus are the functional derivatives 8U[fi]/8fi. 

In frequently employed control algorithms like Krotov 
[13] and GRAPE [HI [22] the time interval in which con- 
trol is exerted is divided into many short sub-intervals 
and the time-dependence of control fields is described 
by step-like pulses that are constant within each sub- 
interval. Consequently, the time-evolution operator U 
becomes a product of propagators for the individual sub- 
intervals that are induced by time-independent Hamil- 
tonians, and the functional derivatives reduce to ordi- 
nary derivatives with respect to the fields' amplitudes in 
the individual sub-intervals. These ordinary derivatives 
can be obtained rather efficiently, which results in the 
astonishing performance of Krotov and GRAPE. This 
decomposition of the control-fields' time-dependence in 
piecewise constant segments, however, often results in 
pulse-shapes that contain high-frequency components. 

In the present letter, an expansion of the field ampli- 
tudes 

fi(t)=J2 a v 9j(t) (1) 
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in periodic functions gj(t) with a given periodicity T is 
employed. Since the period T can be chosen to coincide 
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with the duration over which control is exerted, the as- 
sumption of periodicity does not result in any restriction 
in practice, but due to this periodicity Floquet theory 
[TBI |2"T1 |2"5] can be invoked to determine the propaga- 
tor U. According to Floquet's theorem, the Schrodinger 
equation with a periodically time-dependent Hamiltonian 
with periodicity T has a complete set of solutions 



!**(*)) =e- l£fct |$ fc (t)) 



(2) 



with the quasi-energies £j and periodically time- 
dependent state-vectors \$k(t)) — |$fc(t + T)), which 
are the eigenvalues and corresponding eigenstates of the 
Floquet-operator K, — H — idt- With the solutions given 
in Eq. ^ the propagator reads 



w(*) = S l*fc(*)X**(o)| = 



,—iekt 



<M*)><**(0)| ■ (3) 



In practice, the quasi-energies Ek and vectors |$fc(t)) are 
constructed with the help of the Fourrier-transformed op- 
erator JC = J Q T dt JCe luJt /T of JC. it has a 27r/T-periodic 
spectrum, and the eigenvalues £& and eigenvectors \xk) 
within one Brillouin zone, i.e. interval of width 2tt/T 
determine all solutions 



(4) 



of the Schrodinger equation, where \xkv) is the subvector 
of \xk) associated with the frequency component vuj. 

Since hi is determined completely in terms of the eigen- 
system of K via Eqs. (|3| and Q, the desired derivatives 
of hi with respect to the control parameters are given 
by the corresponding derivatives of £& and \xk)] these, in 
turn, are readily obtained from a perturbative analysis 
of /C. As IC depends linearly on the control parameters 



the linear expansion 



it = ita + y ' K. 



(5) 



with /Co = K\ 
choice of ayy . 



,(°) 



and Kij = d ai IC is exact for any 



The dependence of £& and \xk) on the 
control parameters ay can therefore be described by a 

series expansion with the small parameters — a\j, 
and the expansion coefficients in n— th order are given 
by the n— th order perturbative corrections with /Co as 
unperturbed operator and the perturbations /Cy. The 
explicit form of the derivatives reads 



dsk 

ddij 

d 2 e k 

ddijdcipq 

d 2 \Xk) 
da,ijda pq 



(Xk\ltij\xk) , -r^TT = -ZkTkij\xk) 
—2{Xk\TkijZk7~kpq\Xk)s , 



2 lkTkij1kTk pq \xk) s — lAxk) 



where \xk) is the eigenvector of ACo to the eigenvalue 
£k, the symbol 's' indicates symmetrization with respect 
to {i,j} and {p,q} and (j, = (xk\%ij^lT kpq \xk) s, i-e. 
V = \{Xk\Tkij1lTkpq+Tk P q1lTkij\Xk)- T-k is the pseudo- 
inverse of /Co — £fc ; the inverse restricted to the sub- 
space orthogonal to \xk), and Tkij = K-ij ~ ^frl- 

These derivatives determine the derivatives of the 
propagator hi, and, in turn, of any target functional, i.e. 
property of the final state g(tf) = ht(tf)g(0)W (tf) at 
a given final time tf that is to be optimized. Similar to 
Krotov and GRAPE [TSjrrj] this information can be used 
to improve a pulse sequence as parametrized by the a^, 
and close-to-optimal pulses are obtained in an iterative 
manner. 

The main goal of the present approach was to arrive at 
control pulses with a potentially narrow spectral range. 
Yet, the above expansion of a control sequence also yields 
various additional benefits: 

- A control pulse is parametrized with comparatively 
few parameters; as shown below pulses with less 
than 10 frequency components yield very satisfac- 
tory results, whereas conventional methods with 
piece-wise constant functions typically require the 
optimization of several hundred or thousand pa- 
rameters iiBi rrfi. 



Also higher order derivatives of U with respect to 
the control fields can be obtained exactly, what 
improves convergence as compared to algorithms 
based on first derivatives only [22l [24] . The curva- 
ture is described by quadratically more scalars than 
the gradient, which makes its evaluation intrinsi- 
cally more expensive, but due to the small number 
of control parameters in the present framework, ac- 
cessing the exact curvature is certainly practical. 
If, however, this effort is to be avoided, the (scalar) 
curvature along a specific direction (typically the 
gradient) is also available in a perturbative analy- 
sis with a single perturbation Y]. ■ bijK-ij, where 6^ 
determines the direction along which the curvature 
is to be determined. 

Since the complete time-propagation operator is 
constructed, as opposed to its action on a pre- 
defined initial state, accessing non-linear target 
functionals [25] hardly requires additional effort as 
compared to typical functionals linear like fidelity. 

Boundary conditions, such as smooth switching on 
and off of the pulses, can be realized easily through 
a suitable paramterization in Eq. ([!]). 

Since all dynamics are given by time-dependent ex- 
ponential functions, which can be integrated ana- 
lytically, time-averaged target functionals can also 
be implemented without additional overhead. This 
permits, for example to minimize the deviation 
from a desired dynamics over a finite time-interval. 
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FIG. 1: Logarithmic infidelity and required maximal ampli- 
tude as obtained via 100 steps of iteration for spin ensembles 
with frequency range 2 > uj > —2 (1 > uj > —1) depicted in 
grey (black) and i ma x = 4(7) components in the control pulse 
depicted with *, (x, 0). The inset show the logarithmic 
infidelity achieved with the pulse after 100 iterations as func- 
tion of the resonance frequency of the spin in Eq. d6J. This 
pulse with oi = 0.498 a 2 = -1.904 a 3 = 2.114 a 4 = 0.884 
has been constructed to achieve close-to-unit fidelity in the 
interval 2 > uj > — 2. 



Having established the framework, its efficiency can 
be demonstrated with the paradigmatic example of an 
inhomogeneously broadened spin-ensemble. Each spin is 
described by the Hamiltonian 



7-L 



-A(t)a y with A(t) 



sinf ( ^""^ a,j sin if), (6) 
»=i 



with its natural eigenfrcqucncy uj. All units are taken in 
terms of the natural system frequency so that the scaled 
frequency u>, A(t), a>i and t are unitless. The pulse-shape 
given in terms of the prefactors ai are to be determined 
and the parametrization ensures that the control-pulse 
and its time-derivative vanish at t = and t = tt when 
the control held is switched on and off. The goal is to 
drive a transition from the a z eigenstate |1) to the bal- 
anced superpositions |$) = (|0) + |l))/\/2 of both cr z 
eigenstates. The pulse should achieve this goal for a 
broad range of system eigen- frequencies uj. This can be 
achieved with the help of an ensemble of spins described 
by Eq. ^ with different frequencies uj and the require- 
ment that the pulse realize the desired dynamics for each 
ensemble member. In the present case, an ensemble of 50 
spins with uj drawn from a uniform distribution from the 
interval [—2, 2] was used. Choosing an interval centered 
at implies a transformation to an interaction picture, 
where the average spin is stationary and the part of the 
driving held that oscillates with twice the average reso- 
nance frequency can be neglected |27j . 

Target of control is the averaged fidelity / = (f{(jj)) u , 
where f(uj) — |( , J , Cl j| < l>)| 2 is the fidelity between the fi- 



nal state that is reached with the Hamiltonian H u 
and the average in / is taken over the considered en- 
semble. To ensure the identification of a pulse with low 
intensity, there is an additional penalty function that 
creates an incentive to find weak pulses. Specifically, 
it reads T — J — p J^ f dtA 2 with a non-negative func- 
tion p that decreases as the algorithm reaches its target. 
The inset of Fig. ([!]) depicts the infidelity 1 — f u for a 
pulse with i max — 4 components obtained after 100 iter- 
ations. The logarithmic scale indicates that a fidelity of 
f(oj) ~ 99.999% has been reached in the entire interval 
2 > UJ > -2. 

The lowest order term a% ~ 1/2 corresponds to a pulse 
that achieves the desired dynamics for a spin with uj = 0. 
A pulse with only this component has a maximal ampli- 
tude of 1/2 at t — tf/2. The higher order components 
dj (z > 1) that achieve the robustness of the pulse un- 
der variations of uj necessarily lead to higher amplitudes, 
so that the numerically obtained pulse has a maximum 
amplitude of 2.77 that is attained at t = 0.717 tf. The 
prize for the robustness is thus an increase of maximum 
amplitude by a factor of 5.5. 

Fig. [I] depicts the performance in terms of logarithmic 
infidelity and maximum amplitude of the pulses obtained 
during the iteration. To start with, the pulse is initialized 
randomly with very weak amplitude, and Fig. [T] shows 
how the obtained fidelity increases with growing ampli- 
tude. The symbol 'x' in Fig. [T] correspond to the case 
with spread of resonance frequencies 2 > uj > —2 dis- 
cussed above. The symbol '*' correspond to 1 > uj > — 1; 
not surprisingly, a given fidelity can be reached with 
lower amplitudes than in the former case, and larger fi- 
delities (up to 1 — 10 -7 ) can be reached. On the other 
hand, the performance does not improve if more higher 
frequency components are included in the control pulse 
as depicted with the symbols '0' and '(8)' that correspond 
to pulses with i max — 7 in Eq. ([6]). As it can be seen, 
the amplitude-fidelity characteristic coincides with that 
of the imax = 4 case. That is, there is no improvement 
in going to broader frequency ranges, but only a slight 
decrease in the rate of convergence. 

To gain insight into the dynamics of the spin ensemble, 
it is instructive to inspect three components of the Bloch 
vector, i.e. (l\jU* (t)ajjU (t)|l) (i = x,y,z) with the Pauli 
matrices They are displayed in Fig. ^ for 200 spins 
with resonance frequencies in the interval 5 > uj > — 5. 
The trajectories corresponding the spins in the interval 
2 > uj > —2 for which the pulse had been designed are 
depicted in black; the other trajectories are depicted in 
grey. As it can be seen, during the course of time the tra- 
jectories spread, but there are several instances when the 
individual components of the black trajectories are refo- 
cussed. There are two such events during the propagation 
for the cr^-component and one rather weak focus for the 
o~y component. Only at the end of the propagation at 
t = tf are all the three components refocussed. The grey 
trajectories often deviate strongly from the black ones, 
and, in particular, at tf when the black ones are being 
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FIG. 2: Dynamics of spin components for 500 spins. The gray trajectories correspond to spins with resonance frequencies 
outside the targeted frequency window; the black trajectories correspond to the frequency range within which the desired 
target state is to be achieved. 



strongly refocussed, the grey ones rather spread strongly. 

Whereas the focus here was on the coherent dynam- 
ics of an ensemble of non-interacting spins the present 
approach directly applies also to interacting manybody 
systems, and it can also be extended to dissipative quan- 
tum systems |26j . In particular, in the growing field of 
hybrid quantum systems the present method opens up 
new pathways to control the interaction between compo- 



nents of different types, improve coherence or entangle- 
ment properties, or also, e.g. to improve the accuracy in 
the identification of system properties. 
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